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We consider a weakly-interacting, harmonically-trapped Bose-Einstein condensed gas under rotation 
and investigate the connection between the energies obtained from mean-field calculations and from 
exact diagonalizations in a subspace of degenerate states. From the latter we derive an approximation 
scheme valid in the thermodynamic limit of many particles. Mean-field results are shown to emerge 
as the correct leading-order approximation to exact calculations in the same subspace. 

PACS numbers; 03.75.Fi, 05.30. Jp, 67.40.Db, 67.40.Vs 



Following the observation of Bose-Einstein condensa- 
tion in vapours of alkali atoms , considerable atten- 
tion has been devoted to the behaviour of these systems 
under rotation. Matthews et at ^ reported the obser- 
vation of vortex states in a two-component system, and 
Madison et al. Q demonstrated the existence of vortex 
states in a stirred one-component Bose-Einstein conden- 
sate. Numerous authors have considered the problem 
of rotation in these systems theoretically [5|-[Tl||. Many 
of them focused on the Thomas-Fermi regime of strong 
interactions between the atoms Pj8|jTll| , where the coher- 
ence length is much smaller than the spatial extent of the 
system and where the cloud is expected to exhibit prop- 
erties like those of a bulk superfluid such as ''^He. The 
other interesting limit, to be considered here, is that of 
weak interactions between the atoms, where the coher- 
ence length is much larger than the size of the cloud of 
atoms and where one expects to see properties similar to 
those of superfluid nuclei . 

The problem of rotation in the limit of weak inter- 
actions has also been considered previously. Wilkin et 
al. H] studied the case of effective attraction between 
the atoms and showed that, in the ground state of the 
system, angular momentum is carried by the center of 
mass. Butts and Rokhsar |0] examined the problem of 
effective repulsion between the atoms. They included in- 
teraction effects in mean-field approximation using the 
Gross-Pitaevskii equation and performed numerical cal- 
culations of the moment of inertia of the system. In Ref. 
1^ one of us demonstrated that these systems exhibit two 
additional kinds of condensation associated with the de- 
grees of freedom which come from the multiplicity of the 
degenerate rotating states. Finally, Bertsch and Papen- 
brock jl^ solved the problem exactly within a subspace 
of degenerate states for a limited number of particles and 
angular momenta. 

In the absence of interactions, there is a large degen- 
eracy corresponding to the many ways of distributing a 
given angular momentum among the particles |l^]. The 
challenge is to incorporate interactions in order to lift the 



degeneracy. This can be done either within the mean- 
field approximation or "exactly" by diagonalizing 
the Hamiltonian in an appropriate subspace of the de- 
generate states. In this Letter we establish a connection 
between these two approaches. 

Our starting point is the hamiltonian H = hi + V , 
where 



2M ' 2 



(1) 



includes the kinetic energy of the particles and the po- 
tential energy due to the trapping potential and where 



V 



(2) 



is the two-body interaction between particles, which is 
assumed to be of zero range. Here M is the atomic mass, 
Lu is the frequency of an isotropic trapping potential, and 
Uq — Aiih'^a/M is the strength of the effective two-body 
interaction, with a being the scattering length for atom- 
atom collisions. We assume that a > 0, i.e., we treat 
only the case of repulsive effective interactions between 
the atoms. We further restrict our attention to the limit 
of weak interactions. 



(3) 



where n is the density of the condensed atoms. If the 
system has angular momentum L, this condition allows 
us to work in subspaces of states which are degenerate 
in the absence of interactions. All other states differ by 
an energy of order Tiw, which is much larger than nlJo 
by assumption. Finally, we assume that the cloud of 
atoms rotates about some axis, and that the system is in 
its ground state with respect to this axis, which implies 
that our problem is essentially two-dimensional. In the 
absence of interactions for the harmonic oscillator poten- 
tial in two dimensions the single particle energies are 



e = [2nr 



(4) 
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where Ur is the radial quantum number, and m is the 
quantum number corresponding to the angular momen- 
tum. In the lowest energy state of the system all the 
atoms are in states with rir — 0, and with m being zero 
or having the same sign as the total angular momentum. 

Let us describe briefly the two methods which have 
been used to attack the problem of rotation, starting with 
the mean- field approach Within this approxima- 

tion the many-body wavefunction with N particles and 
L units of angular momentum ^'ijv(ri, r2, . . . , rjv) is a 
Fock state simply expressed as a product of single parti- 
cle states: 

*L,7v(ri, r2, . . . , Tat) = ^'(ri) x ^(ra) . . . *(rjv). (5) 

It is natural to expand the single-particle states ^'(rj) in 
the basis of the harmonic-oscillator eigenstates ^m{^i) 
with angular momentum mh along the axis of rotation: 

oo 

\E'(r,) = ^ c,„$„(ri), (6) 

m=0 

where Cm are variational parameters. The quantity |cmp 
gives the probability for the occupancy of state $m . One 
imposes two constraints on the parameters Cm- The 
wavefunction must be normalized, 

m 

and the expectation value of the angular momentum per 
particle must be fixed, 

J2Mcm\^ =L/N. (8) 

m 

The expectation value of the interaction energy V in the 
state given by Eq. (||) is 

{V) = ^N{N-l)UoJ \^\Ur. (9) 

Minimizing this energy with respect to the Cm subject to 
the constraints of Eqs. (0) and (||), one obtains the mean 
field energy and corresponding wavefunction. 

Turning to the "exact" approach to this problem |^,|lO| , 
it is convenient to work in the basis | A^o, -^ij ^2, ■••) where 
is the number of particles with angular momentum 
hm. The constraints on particle number and angular 
momentum now take the form: 

J2 ^N;J2 = (10) 

rn m 

It is convenient to write the interaction energy as 

V = ic/o X! ^^'0'''''' aja^jUkai, (11) 

where ak and aj. are annihilation and creation operators 
respectively, and 



I.J.kJ = J ■I':(r)<I'-(r)*i(r)*,(r)'ir 

when i + j — k + I and zero otherwise. The final step 
is to diagonalize the matrix V in the space of the states 
|iVo, A^i, iV2, ...) which satisfy the constraints of Eq. ([l0|). 

Before turning to a comparison of these methods, it 
is useful to consider their advantages and disadvantages. 
On the one hand, the mean-field approximation gives the 
energy in the asymptotic limit L and N oo with L/N 
fixed. This limit is not directly accessible to the exact 
method since the dimensionality of the space increases 
dramatically with increasing L and N. On the other 
hand, the exact method provides the entire spectrum of 
excited states and a lower ground state energy than the 
mean- field result. The mean- field assumption of a simple 
product wavefunction ignores correlations between the 
particles which are contained in the exact diagonaliza- 
tion. The mean-field energy offers only an upper bound 
for the true ground state energy. To be more specific, the 
leading term in the ground state energy of the system 
must grow like N'^ for the interaction considered. As we 
show below, the energy obtained by the exact method in 
the asymptotic limit is bounded from below by an expres- 
sion which has the same iV^ coefficient as the mean-field 
result. Thus, the exact energy differs from the mean-field 
value by terms of order N or less. In addition, we will 
suggest simple methods for computation of the order N 
correction to the ground state energy in the asymptotic 
limit. 

The spirit of the mean field approximation is that the 
wavefunction is composed of a large number of states 
with occupation numbers in the vicinity of iVjcmP but 
that these states are of measure zero in the full Hilbert 
space considered. Up to a staggering related to the 
phases of the c^, the wavefunction should then be a 
smooth and differentiable function of occupation num- 
ber for single particle states containing a finite fraction 
of the particles. Such states are amenable to simple treat- 
ment. States with Cm = in a mean field description will 
require special treatment. Consider the particular case 
L — 2N using a Hilbert space < m < 4. In mean field 
theory it is found that the only non-zero coefficients are 
Co, C2, and C4. Initially retaining only these states, we 
will consider the basis of states 

= |fc,0,iV- 2fc,0,fc). (13) 

Write the wavefunction as 

\L = 2N,N)=^{-l)'^ijk\k), (14) 

k 

where the fc-dependent phase has been chosen to mini- 
mize the resulting energy. The eigenvalue equation has 
the form 

Vk,k1pk - Vk,k-l^k-l - Vk^k+li^k+l = S2N,N'tpk, (15) 
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where Vk,k' = (^1^1^') is the matrix element of the inter- 
action between the states and \k') given by Eq. (13). 
The two-body interaction V connects only states for 
which k — k' = ±1, or 0. 

Assuming that ipk is smooth and differentiable over a 
suitable range, we expand the wavefunction as 



V'fcii = -0?= ± dktpk + ^dlil^k (16) 
to obtain the eigenvalue equation 

- ^(Vfe,fc-i + Vk^k+i)dl^pk + Vcff^pk = S2N,N'ipk, (17) 
where the effective potential is 

VMk) = Vk,k - Vk,k-i - Kk+i- (18) 

The leading term in the energy is of order N'^ and is 
simply the minimum value of Vcf[{k). This minimum 
occurs at 



nn = 



\/d{VkM-i+Vk,k+i) 



0.1192iVwo- 



(23) 



N 



16V6- 28 
64\/6- 109' 



(19) 



The occupancy of the m = and m — 4 states is ko/N 
in the asymptotic limit; the occupancy of m = 2 is 
1 — 2ko/N. This is precisely the result found in Ref. 

within the mean-field approximation. Numerical ex- 
ploration for TV > 10^ indicates that the energy of the 
minimum of the effective potential, Ves{ko), can be de- 
scribed by the approximate form 

^2^,^ ~ (0.1772568957V2 - 0.366043iV -I- 0{N°)) vq, 

(20) 

where vq — Uq / |<i>o|'* dr, with $o(r) being the ground 
state of the harmonic oscillator. As expected, the N'^ co- 
efficient in this expression is identical to the correspond- 
ing term in the mean field |^^. In general, the properties 
of the minimum of Vcs{k) reproduce the leading N'^ term 
of mean field theory. 

Expanding Vcs{k) around fco, we find that the effective 
potential has the form 



Vcs{k) = i2N,N + -d{k- ko^vo, 



(21) 



where £2N,n is given by Eq. ( poD and where the constant 
d « 0.3732 is necessarily positive and of 0{N'^). The 
leading, "kinetic energy" term in Eq. ( pT| ) is positive with 
Vk,k-i + Vk,k+i ~ 0.0381iV2. That this term is of 0{N'^) 
is a consequence of the fact that y is a two-body opera- 
tor. The ground-state energy for this harmonic oscillator 
problem is readily found to be 



£2N,N = £2N,N + 



(22) 



with 



The effective oscillator length corresponding to Eq. ( |l7| ) 
is proportional to {\Vk,k-i + Vk,k+i\/dY/^ and is of order 
N^^^. This demonstrates the consistency of our assump- 
tions: The ground state wavefunction ipk is a Gaussian 
centered at fco with a width of order N-^^^. It contains 
significant contributions from a large number of states, 
which, however, is only a vanishingly small fraction of 
all the states in our Hilbert space. Combining Eqs. ([2^), 
(p^), and (^) we obtain a ground state energy of 

£2N,N = (0.177256895A^2 _ o.306406A^ + 0{N°)) vq. 

(24) 

The general conclusion from this analysis is that the 
correct 0{N'^) contribution to the energy can be obtained 
by determining the classical minimum of the effective po- 
tential energy, Vcs_{k). Quantum fluctuations lead to cor- 
rections of order TV. Since the kinetic energy term is pos- 
itive, these corrections are necessarily positive and this 
leading term thus provides a lower bound on the ground 
state energy. By contrast, the mean-field approach leads 
to minimizing the same classical potential but provides 
an upper bound on the energy as a consequence of its 
variational character. Thus, the ground state energy is 
trapped between upper and lower bounds with the same 
leading term. Thus, we conclude that the 0{N'^) contri- 
bution to the ground state energy obtained by mean field 
theory is exact in the large N limit. 

While the space considered in the above example was 
severely limited, this method can be generalized to in- 
clude as many states as necessary. Our analysis holds 
provided that the states included are occupied by a finite 
fraction of the total number of particles. If this is true, 
roughly A^^/^ states for each of the variables in the multi- 
dimensional space will make significant contributions to 
the wave function and the approximation of differentia- 
bility is valid. One can then derive an eigenvalue equa- 
tion of the form of Eq. (|l^) , of multi-dimensional coupled 
harmonic oscillators. A crucial assumption in establish- 
ing the asymptotic validity of the mean field result is the 
positivity of the "inertial" parameters such as the coef- 
ficient of d^^k in the case of Eq. (p^). We have verified 
that this is true in a number of specific cases, and we 
believe that it is true in general. 

Let us now return to the states m = 1 and m = 3 in 
our example. These states are not macroscopically oc- 
cupied, i.e., the corresponding ci and C3 are zero in a 
mean field calculation. The expansion of Eq. (^6|) can- 
not be used for their description, and their incorpora- 
tion requires some care. To include their effects, we ap- 
proximate \L = 2N,N)q by the single state \k = ko) of 
Eq. dl^). Starting with this state, we build a basis of 
states containing ni particles with m = 1 state and 
particles with m = 3, being careful to conserve parti- 
cle number and angular momentum. The effect of these 
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states with odd m values can then be obtained by diago- 
naUzing V in this basis. If < (ni +713) < 2Q, there will 
be (Q+l)^ states in this basis. Convergence of this calcu- 
lation to arbitrary accuracy can be obtained by choosing 
Q to be large but none the less of order N'^. This con- 
firms the fact that the states m = 1 and m = 3 arc not 
macroscopically occupied. In the present case, the lowest 
eigenvalue of this matrix converges rapidly to an asymp- 
totic value of — 0.1812A^wo. The linear contribution to 
the energy from such states is general. An improved de- 
scription of |L = 2N,N)q would lead to changes in this 
energy of order iV°. 

We are now able to obtain the asymptotic ground state 
energy for the case L = 2N correct to order N by com- 
bining this correction for the odd m states with the result 
of Eq. (pi. We find 



-2N.N 



(0.177256895iV2 - 0.487576iV + 0{N°)) vq. (25) 



This result is in remarkable agreement with the expres- 
sion 

£2N,N ~ (0.177256752iV2 - 0.487848iV + 0(iV")) vq (26) 

obtained by direct numerical diagonalization in the same 
truncated space < to < 4 for values of 100 < N < 220. 
Note that the case N = 220 requires diagonalization of 
a matrix of dimension 305990. The small discrepancies 
obtained are a consequence of the limited range of N for 
which exact diagonalization is possible and the largest 
values of Q considered in our approximate calculation. 
These restrictions are absent in calculations within the 
same Hilbert space < to < 4 for L = N/2. It is then 
possible to consider N as large as 600. All states have 
an occupation number of order N, and the final matrix 
diagonalization in our approximate treatment is not re- 
quired. In this case, the quadratic terms in the ground 
state energy agree to 10 significant figures, and the linear 
terms differ by 1 part in 10^. 

Additional evidence comes from exact diagonalizations 
(with no truncation in m) for the case L = 2N for 4 < 
N < 24. The resulting energies are well described by 



£2N,N ~ (0.17497V2 - 0.5206iV + 0.0028)wo. 



(27) 



As argued above, mean field theory suggests that only 
states with even to are macroscopically occupied in the 
ground state for L — 2N. This is supported by exact 
diagonalizations for the same range of N from which we 
have excluded single-particle states with odd m. The 
energy is then found to be 



2N.N 



(0.1755Ar2 _ o.2980iV - 0.0905)wo. 



(28) 



Comparing the coefficients of and N in Eqs. ( |27| ) and 
( p8| ) , we see that only the latter is materially affected by 
this truncation. This is consistent with our assertion that 
states which are not occupied macroscopically contribute 
to the energy at order N. 



In summary we have adopted a "correspondence prin- 
ciple approach" to describe repulsive two-body interac- 
tions between the particles in a Bose- Einstein condensate 
by assuming that the wavefunction in the N 00 limit 
is a smooth and differentiable function of the occupation 
number. This approach leads us to the result that mean 
field theory provides an exact determination of the iV^ 
term in the ground state energy and permits determina- 
tion of the contribution to this energy of order N. Butts 
and Rokhsar have arrived at a similar conclusion . Ev- 
idently, the same arguments can be used to describe the 
spectrum of excited states. The present methods would 
appear to be a useful and practical supplement to mean 
field methods whenever these are of value. 
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